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THE ABSENCE OF RAPID X-RAY VARIABILITY IN ACTIVE GALAXIES^ 

Allyn F. Tennant^ and R.F. Mushotzky 

Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 


We have searched for variations on time scales ranging from minutes to 
several hours In the X-ray flux from 54 observations of 38 active galaxies. 

Our sample Is composed mostly of Seyfert I galaxies but also Includes radio 
galaxies, NELG's, BL Lacs and 3C 273. Only NGC 6814 varied on time scales as 
short as 100 sec. No other source was observed to vary with a time scale of 
less than 12 hours. We conclude that large-amplitude short-term variations 
are not a characteristic of the X-ray emission from active galaxies. Upper 
limits on Oj/I ranged from Z% for Cen A, 5% for NGC 4151, to 20% for sources 
giving 1 ct/sec In our detector. Three objects NGC 3227, NGC 4151 and MCG 5- 
23-16 show variability consistent with a time scale of ' 1 day. 

We consider ways to reconcile the rapid variability seen for NGC 6814 
(and NGC 4051) with the general stability observed for the other objects. If 
the lack of variability Is due to the growth of a new, constant source of soft 
photons (such as thermal re-rad1at1on from dust) then the X-ray plasma must be 
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much larger than the size Implied by the variations seen In NGC 6814. Other 
possibilities are that the "originally lumpy" accretion flow stabllzes or that 
NGC 6814 Is different from other galaxies. 

Subject headings — Galaxies; nuclei - Galaxies: Seyfert - X-ray: sources 
Iprom a dissertation to be submitted to Graduate School, University of 
Maryland, by Allyn F. Tennant In partial fulfillment of the requirements for 

the Ph.D. degree In Physics. 
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I. INTRODUCTION 

It is currently popular to assume that the X-rays from active galactic 
nuclei a*e produced deep in the potential well near a super massive object 
(see Rees, Begelman and Blandford 1981 for recent discussion). If the X-ray 
flux from this dense gravitationally confined plasma is observed to vary then 
the shortest time scale for variability will be on the order of the light 
travel time across the most stable orbit in a Schwarzschild geometry. This is 
given by 

At - R/c ' 10 GM/c^ - 50 sec Mg 

where Mg Is the mass of the central object in units of 10® Mq. Since Mg is 
expected to be in the range of 1 to a 1000 the relevant time scales range 
from ' 1 minute to ' 1/2 day. Lack of variability on these time scales could 
be an indication that the X-ray emitting plasma is not gravitationally 
confined. 

Recently, Guilbert, Ross and Fabian (1982) have discussed a model which 
accounts for the simple power law energy spectrum observed from these 
sources. In their model a region of hot electrons produces X-rays and cools 
via inverse Comptcnization. Since the cooling time is much shorter than the 
period in which a good spectrum can be obtained the observed spectrin will be 
a time average. Sources powered with this mechanism should show a low 
amplitude flicker with a time scale of minutes. 

On the observational side, before the launch of the HEAO observatories 
only two active galaxies, Cen A and NGC 4151, were bright enough to have had 
their X-ray flux accurately measured in less than 12 hours. Both show 
vai iv'.hil ity on time scales less than one day. Winkler and White (1975) ';-aw a 
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factor of 1.6 change during a 6 dey observation of Cen A. The angular size 
Implied by such a change Is less than 1 marc sec. De1va111e» Eptein and 
Schnopper (1978) found evidence of a 25% Increase In Cen A during a fev> hours* 
confirming rapid variability. Lawrence* Pye and Elvis (1977) have shown that 
Cen A Is continuously variable In Its X-r^ flux over extended periods of 
time. 

The second brightest c- tragalactic X-ray source NGC 4151 also has a long 
history of observations of X-ray variability. Mushotzky* Holt and Serlemltsos 
(1978) reported a flare with ft factor of 2 change ob'erved over 1.5 days. 
Lawrence (1980) using Ariel 5 data extending over several years finds that 
such flare-like events are rather common. Although the absorption column does 
vary In NGC 4151 and causes some of the observed flux variations (Barr et a1. 
1977) there are no Indications of changes In the absorption or spectral Index 
during the more rapid changes (Mushotzky, Holt* Serlemltsos 1978; Baity et a1 . 
1982). 

Since the two brightest active galaxies are highly variable on time 
scales as short as one day. It was expected that this should be a property of 
all active galaxies. In the pre-HEAO days It was difficult to detect rapid 
variability due to small collecting areas and low on source duty cycles. 

Ariel 5 results (Pounds 1979; Marshall, Warwick and Pounds 1981) support the 
possibility that all active galaxies (occasionally at least) show variability 
on time scales as short as one d«y. However, It was Impossible for Ariel 5 to 
measure time scales shorter than one day for weak sources. 
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In this paper we report on the results of a search for rapid variability 
from active galaxies using the HEA0>1 A2* experiment (Rothschild et a1. 

*The A2 experiment on HEAO-1 Is a collaborative effort led by E. Boldt of 
and G. Garraire of CIT, with collaborators at 6SFC, CIT, JPL and UCB. 

1979). Since we are Interested In the most rapid time scales we have used the 
pointed data. In this mode the spacecraft observed the source of Interest more 
or less continuously for a period lasting from 3 to 18 hours. There were 54 
pointings which gave data useful for this analysis. The high on source duty 
cycle combined with the low background, large bandwidth and large collecting 
area resulted In high sensitivity to variations with time scales ranging 
from ' 1 minute up to ^ 3 fours. Longer time scales will be discussed by 
Mushotzky et a1. 1982. 

II. METHOD OF ANALYSIS 

For this analysis we used data from one Xenon detector, HED 3, which 
means that we are most sensitive to X-rays from 3 to 15 keV. We selected data 
with low veto rates (electron quiet) and low Mcllwaln L values (see below). 

We also excluded times when the earth or earth's atmosphere was In the field 
of view. These data selection criteria resulted In ' 70% of the data being 
rejected for further analysis. We subtracted a constant background rate and 
then corrected for effective exposure to the source to obtain the source 
flux. Errors In the source flux were assumed to be purely statistical with an 
additional term due to the uncertainty In the absolute level of the 
background. 

2 

It was then possible to compute x testing the assumption that the source 
did not vary during a given observation. Since the probability of observing 
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any given value of reduced (that Is pe>^ degree of freedom denoted by xj) 
Is a strong function of the number of degrees of freedom N, we computed 


X • (xj - D//OT. 


( 1 ) 


For large values of N one can convert X Into a probability of X being less 
than or equal to this value by chance, using Table 26.1 of Abramowltz and 
Stegun (1970). We have found that for N ranging from a few to Infinity using 
X » 1.3 gives a 90% confidence upper limit for detection of variability. The 
usefulness of X Is lost If ^ thus for X values of less than 1.3 we 
define a 90% confidence upper limit by computing how much the source would 
have had to vary to generate an X value of 1.3 exactly. 

We have separated the variance Oj In the source intensity I from the 
variance Introduced by counting statistics using 


Ot - 

= ( 4 ) 


"Xn 


- 1 ) 


H/Z 




( 2 ) 


In equation (2) Oj Is the error In Ij for the data bln {see Appendix A); 
where 0 Is a measure of the source variability and is a function of the time 
scale over which one averages and the "true" timescale for source 
variability. Since Cj Is normalized by dividing by the source Intensity the 
resulting number Is dimensionless and so measures an Intrinsic property of the 
source (Independent of distance). We have numerically tested equation (2) by 
generating "phantom" data with a -known amount Df "source" variance and adding 
In the correct amount of counting statistics variance. We have found that for 
X > 1.3 equat1r*n (2) did allow us to separate "source" variance from photon 


noise. 
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For a variable source the value of $ Is a function of the size of the 
time bins one used to calculate x • example If the bln size is shorter 
than all relative time scales a will be at a high and roughly constant 
value. However, for longer bln sizes one starts to average more than one 
"event" Into a single bln and the value of a will decrease. For a plot 
of a vs bln size for Cyg X-1 see Dol (1980). An early version of this 
analysis was done by Boldt et al. (1975). 

It will be useful below to consider the value of a for the case where the 
Intensity Is uniformly and randomly distributed between a low of Iq and a high 
of Ij. When the bln size Is smaller than the transition time, we find 

2 h ■ ^0 Al . 

^ 717 I, + ■ <r>717 

where Al * and <I> Is the average. For bln sizes equal to or longer 

than the variation time scale, one starts to average over fluctuations which 
will cause the observed value of a to be reduced from the value given above. 
Thus on a a vs. bln size plot, the location at which the curve rolls over Is 
related to the time scale of the dominant source of variance In the data. 

III. SOURCES OF ERROR 

In any search for variability one expects to find a few events that are 
statistically quite rare but nonetheless are not "real" events. This can 
arise from two quite different causes. First, extremely rare non-X-ray 
related events can occur in the detector giving rise to spurious 
variability. Secondly, there can be a small amount of residual noise In the 
detector which gives rise to an excess variance (above Poisson noise). The 
residual noise can be totally invisible in a single observation but can, in 
the case of a large number of observations, give rise to more than the 
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expected number of "3o" detections of excess variance. In this section we 
consider both types of errors. 

A. Non-X-Ray Events 

Twice out of a total of 102 observations, once observing NGC 526a and 
once observing NGC 4593, we observed a short duration "flare". For the case 
of NGC 4593 the flare Increased the "source Intensity" by greater than a 
factor of 5 on a time scale of 500 sec. In both cases the counting rate In 
the detector not observing the source {see Tennant et a1. 1981, Appendix A) 
also Increased. Therefore, In neither case can the variability be considered 
to be coming from the source. We believe these events occur during times when 
the particle flux has an unusual energy spectrum and a flux low enough not to 
trigger any of our flags indicating particle events. 

B. Aspect 

The standard aspect solution provided by Goddard can be In error by up to 
.15° during pointed observations. We ob^-alned ^ 25 aspect solutions from the 
A-3 experimenters who generate better solutions. We f(und that for most cases 
the Goddard supplied aspect was accurate to better than .05°. An error In 
aspect of this magnitude will produce a o^/I of less than 2 % for our 3° x 3° 
field of view. Since we have used the A3 solution for most of the brighter 
sources and the more variable weaker ones we believe that the residual 
variations found for some of these sources Is not related to aspect 
uncertainties. 

C. Mcllwaln's L Parameter 

The background counting rate In our detector Is related to the value of 
Mcllwaln's L parameter (Mcllwaln 1961). Roughly speaking Mcllwaln L Is the 
distance from the center of the Earth of a trapped particle when It crosses 
the geomagnetic equator. Observations during points at "blank sky" (that Is 
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times when no source of flux > 3 x 10*1^ ergs/cm^sec was in our field of view) 
showed that the excess counting rate aC, which is correlated with L,, fits the 
model aC - 2(L-1) ct/sec for 800 cm^ of detector. This linear fit is accurate 
to 20% for L ranging from 1 to 1.5 (units of earth radii). Since HEAO-1 was 
in low earth orbit, L was generally below 1.2 but occasionally values of 1.5 
or more were encountered. We decided to discard all data with L values of 1.2 
or greater. 

This has become one of the standard data selection criteria which was 
used to analyze all pointed observations of active galaxies and blank sky. 

This flag caused >• 30% of previously good data to be discarded but decreased 
the excess variance by a factor of greater than 60%. If we assume that the 
remaining L values are uniformly distributed bet^^en 1 and 1.2 we can then 
apply equation 3 to estimate that the excess variance due to the residual 
Mcllwain L effect is roughly 

0^^ - AC/(/12 A) 

2 

= .20 ct/cm-ksec 

where we have assumed that the mean area A during pointed observations Is 
close to 600 cm2, we note that AC//12 = .12 ct/sec Is only 0.8% of the total 
(internal and diffuse X-ray) background counting rate of ^ 14 cts/sec. 

D. Background Fluctuations and Source Confusion 

The background Is uncertain to about 3% at any point on the sky In our 3 
x 30 beam. This is due to weak unresolved sources below our source confusion 
threshold. This absolute uncertainty In the background flux Introduces two 
types of errors Into our results. First, this "sky confusion" noise, also 
called fluctuations In the diffuse X-ray background, introduces an Irreducible 
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uncertainty in our absolute flux level of - 0.5 ct/cm^-ksec. For this reason 
Me have not considered any sources with fluxes less than 1 ct/cin^>ksec. 

The second way background fluctuations can affect our results has its 
origin in the fact that HEAO-1 was not always pointed at the source of 
interest during a given observation but rather drifted over a band that 
allowed offsets up to 1/2®. Although the counting rate has been corrected Tor 
mean exposure, we find that the different sky backgrounds observed by the 
detector as it looks at slightly different parts of the sky can cause an 
excess variance to appear in the data. To estimate the amount we first note 
that if the detector explores the entire area of the sky that it is allowed 
to, then only 50% of the total solid angle viewed during the observation is in 
the field of view at all times. We expect the variance for overlapping fields 
to be less than the ' 3% observed for independent 3® x 3® fields. To first 
order, background fluctuations will introduce uncertainties of up to 0.32 
ct/cm^-ksec (1.5% of total detector background). Therefore, including both 
residual particle events and background fluctuations, the maximum excess 
variance we expect Is ^ 0.4 ct/cm^-ksec. We note that the observed excess 
variance can be much smaller than this If the sky background flux Is 
relatively uniform In the observed region or If the detector does not explore 
the full region of sky available to It or if the value of Mcllwaln's L 
parameter was unusually constant for a given observation. 

IV. RESULTS 

A. Variance 

In Table 1 we have surmnarized the results of our tests for source 
variability (the list Is alphabetic by source name). Columns 7 and 8 give the 
value of Oy/I for two different bin sizes. Also Included In Table 1 is the 
day the observation was made where day 1 = Julian day 2441344.5 = Jan 1, 1977, 
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and the mean flux <I> for the observation. Here and throughout the paper we 
report fluxes with units of ct/cm^-ksec. A flux of 1 ct/cm^-ksec In the Xenon 
detector corresponds to 1 x 10*11 ergs/cm^-sec In the 2-10 keV band or 1.8 x 
10*11 ergs/cm^-ssc In the 2-20 keY band assuming a power law spectrtm with 
energy slope a ■ 0.7. Again we point out that <I> Is uncertain by 0.5 
ct/cm^-ksec due to the 3t uncertainty In the sky background. 

We have plotted e, obtained using a 328 sec bln size vs. source Intensity 
In Figure 1. Since $ Is computed using the entire observation. Figure 1 Is 
relevant for time scales ranging from 300 sec to about 3 hours. The solid 
line In Figure 1 corresponds to a constant standard deviation of 0.4 
ct/cm^-ksec which was estimated In the last section to be the upper limit for 
systematic errors. Thus any positive detections below this line could be due 
to residual sky plus detector noise. Orly three objects, NGC 6814, MGC 4151 
(one time) and Cen A (one time), lie significantly above the curve. 

We have examined a down to a bln size of 5.12 sec (a detector readout 

time), uue to counting statistics nur upper limits at 5.12 sec are about two 

times larger than the corresponding values at 328 sec. Since there are fewer 

positive detections of variability at 5.12 sec than at 328 sec we conclude 

that there Is not a large source of variance with time scales of less than 328 

sec. We conclude that large- amplitude, o^/I > lOt , short term variations on 

4 

timescales 5 < t < 10 sec are not a characteristic of the X-ray emission from 
active galaxies. 

As one goes to longer bln >l;es then most of the residual non-source 
related noise Is averaged Into a single bln. Since the spacecraft look 
position generally samples the entire area available to It during one 
spacecraft orbit, sky noise Is greatly reduced for a 90 min bln size. For a 
90 min bln size one also averages over the particle background for one 
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Spacecraft orbit. A plot of S vs. source flux for a 90 min bln size Is 
sho*<n In Figure 2. 

We have attempted to determine our systematics for the 90 min bln size by 

performing the equivalent test for variability during observations of "blank 

sky". Of 49 observations of "blank sky", 17 were not constant at the 901 

confidence level. Since for all active galaxies, except Cen A, the background 

gives more counts than the source, we would expect a similar fraction of the 

active galaxies to show variability. Thus we expect 18 t 4 positive 

detections when In fact the observed value was 21. Thus, to first order, the 

majority of our positive detections only Indicate that there Is a small amount 

2 

of excess variance, beyond counting statistics, not Included in our x test 
for variability. 

If the same excess variance accounts for both the Increases number of 
positive detections for the "blank sky" and the active galaxies then the 
magnitude of the variance should help us decide If any detected variability Is 
real or not. In Figure 3b we have plotted a histogram of the standard 
deviation Og In the background counting rate B normalized by dividi by B for 
a ' 90 min bln size. The white areas represent upper limits and the shaded 
boxes correspond to observed values. We have computed the excess variance for 
the entire set of observations by adding together the excess variance from 
each cbservatlon. For the 90 min bln size this gave o„/B of 0.71 which Is 

D 

what one would expect If most of the excess noise was due to background 
fluctuations only. Thus the observed distribution is consistent with the 
expected excess variance being due to "confusion noise". 

For r . •♦Ions of active galaxies we expect a slightly different 
distribw. riist, when the detector is pointed at an active galaxy the 
counting rate is higher (by definition). This means that our sensitivity to 
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non-source related variance is smaller during observations of active 
galaxies. This Is shown by the fact that the distribution of upper limits has 
moved to the right (Figure 3a). Since the flux from active galaxies was first 
area corrected for effective exposure, to convert to Oj/B we 
multiplied 0^1 by S/B where S Is the source counting rate (source flux times 
mean area). This means we have effectively Ignored corrections due to 
changing area as the detector passes over the source. 

Based on the results of Figure 3, we divide the sources Into 3 groups: 

1) Those with o^/B >2%, 2) Those with 1.5% < o./B < 2% and 3) those with o./B 
< 1.5%. Sky fluctuations cannot account for the variability seen for the 3 
sources In group 1, whereas any variability seen In group 3 can be totally due 
to sky noise. Sources In group 2 are suspect and need further checks. In the 
next section we will consider the sources In groups 1 and 2 1n greater detail 
B. Light Curves 

The most variable source was NGC 6814 (Figure 4). As was discussed 
before In Tennant et al. (1981) this source Is highly variable and shows 
factor of two changes In flux on all time scales down to a few minutes. Table 
1 shows 8 maintains a value near 43% for all the bin sizes considered up to 90 
min. Using the published light curve from the Einstein Observatory Imaging 
Proportional Counter (IPC) observation of NGC 6814 (Tananbaum 1980) we find 
that 8 ' 21% for a 3 hour bln size. This could Indicate that either the 
dominant source of variance has a time scale of hours or that 8 decreased 
perhaps due to a decrease In flux between the HEAO-1 A2 and the IPC 
observations. 

In Figure 5 we Illustrate the difficulties In determining the time scale 
for a poorly sampled observation. The light curve Is for our most variable 
observation of NGC 4151. In the top panel we fit all 6 points to a straight 
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2 

line and calculate a x oF 11.12 for 4 degrees of freedom (dof) which Is an 

p 

unacceptable fit to the data. (This compares to x ■ 39.80 for 5 dof 

calculated assuming the source was constant). In the second panel we have 

2 

tried a model of a step function. For this ca»ie x “ 5.39 for 4 dof which Is 

acceptable. Although It appears (to the eye at least) that the transition 

occurred over the 3 bir>$ nearest where we located the step we point out that a 

small systematic deviation lasting 45 min could easily be Introduced by a 

spacecraft orbital effect. In the 3rd panel we have again fitted a straight 

line to the data but In this case have thrown out the first point. For the 
2 

third case x 4,76 for 3 dof. If the source continued to brighten at the 
rate Indicated In the 3rd panel then It would double Its Intensity In 12 
hours. Thus statistically we cannot distinguish between a 12% flux Increase 
on a time scale of 15 min or the start of a 12 hour flare. We prefer the 
latter Interpretation since It Is consistent with previous observations of 
variability from NGC 4151 (Mushotzky et al. 1978; Lawrence 1980). This 
ambiguity as to the correct model Is not resolved as one goes to shorter bln 

sizes, for this case at least, since the smaller bln size generates no new 

Information. 

The third most variable source Is Cen A (Fig. 4). We note that HEAO-1 
scanning data for Cen A show It to be slightly extended, presumably due to 
weak nearby sources. Marshall and Clark (1981) have reported a nearby source 

which will appear In our field of view. This source confusion can explain 

most of the observed variability. Notice that there Is no Indication of a 
linear charge over the 12 hour observation. 

We now consider the sources with only weak evidence for Intrinsic 
changes. First consider the light curve for H1649-595 In Figure 4 (H1649-595 
was tentatively Identified with NGC 6221 by Marshall et al . 1979). This light 
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curve clearly shows that most of the excess variance Is due to short duration 
“flares" near the end of the observation. For this source the three high 
points are due to i bright confusing source. This was clear when we noticed 
the source “flared" only when the detector was looking at a certain region of 
the sky. Since all light curves were looked at by eye. we know that no other 
source showed such flare like events. We have Included H1649-595 In our 
sample as a warning that “statistically real" variability could have many 
causes. 

Both NGC 3227 and MCG 5-23-16 showed evidence for a linear Increase In 
flux. If such a rate Increase continued then the source f'ux would double In 
10 hours for NGC 3227 and 28 hours for MCG 5-23-16. Although the "time 
required to double the source flux" Is not a measure of source "size" It Is 
still a very useful time scale. For example, it can be applied to the Fabian 
and Rees relation (1979). 

Finally, we come to the MCG 8-11-11 observation shown in the center panel 
of Figure 4. In this figure we see small peaks near 2.1, 3.9 and 5.5 hours. 
Since these points are ^ 1.5 hour apart, i.e. one spacecraft orbit, we suspect 
that these are not real events. A Oj/B of - 1.6% confirms that the 
variability seen can be due to confusion noise. If one ignores the short term 
variability then a weak linear increase is seen which will double the source 
flux in a few days. We are unable to judge the reality of such a trend. 

It is interesting to note that all the variable objects in our sample, 

NGC 3227, NGC 4151, NGC 6814 and MCG 5-23-16 are low luminosity objects with 
< 3 X 10^3 erg/sec. This agrees with the HEAO-1 longer time scale data 
(Mushotzky et al . 1982) which show that lower luminosity sources have a 
greater probability of being variable. In addition, with the exception of 
NGC 6814, the observed variability is consistant with a time scale of about 
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one dey. This result agrees with the Ariel 5 observation (Marshall, Warwick 
and Pounds 1981) that many sources show a one day time scale. 

Variability on longer time scales will be reported by Mushotzky et al. 
(1982). However, we can compare the count rate for different pointed 
observations of the same source to obtain some information about longer time 
scales. We find that Cen A and NGC 4151 show large changes over 6 months. MK 
509, observed four times in two weeks, shows a large amplitude charge. This 
was reported by Dower et al. (1980) using the HEAO-1 A3 data. Our data for 
days 657 to 666 show that the time required for aL = <L> is 22 days. Thus we 
confirm the general nature of the variability reported by Dower et al. but 
indicate a slightly longer time scale. 

V. DISCUSSION 

Our observations lead us to conclude that the active galaxies in our 
sample do not vary on time scales ranging from minutes up to a few hours 
(upper limits from Table I) with Oj/I > .1. In our sample only NGC 6814 
showed evidence of variability on time scales of less than 10^ sec. NGC 4151, 
the source with the second largest value of Oj, had a b of less than 
one-eighth the value for NGC 6814. To our knowledge the only other active 
galaxy which clearly varies on a time scale of less than 10^ sec is NGC 4051 
(Marshall et al. 1981). In this section we try to reconcile the rapid 
variability observed from NGC 6814 and NGC 4051 with the lack of similar 
variability from other galaxies. Although we will consider ways that NGC 6814 
is similar to other objects, we note that it could be quite different. 

A. A New Class of Active Galaxies? 

In this section we will consider NGC 6814 and NGC 4051 together and 
search for any characteristic that sets them apart from other active 
galaxies. Rieke (1978) found both of these objects had low, but not unusually 
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1ow» IR luminosities. When he compared his measurements with those from Stein 
and Weedman (1976) the two objects that differed the most were NGC 4051 and 
NGC 6814, although Rieke points out that Stein and Weedbnan's value for 
6814 was probably In error. The difference may be considered as weak evidence 
for IR variability. More recently. Glass (1979), reporting on IR observations 
of active southern galaxies, notes that from an IR point of view NGC 6814 (and 
NGC 3783?) are only marginally Seyfert like. 

Concerning optical correlations Yee (1980) reports on the very strong 
correlation between the luminosity In Hs and the "nonthermal" luminosities for 
quasars and both broad and narrow line active galaxies. NGC 6814, NGC 4051 
and NGC 3227 lie near each other on the correlation but are a factor of 50 
weaker than the next strongest Seyfert I gala)^ In Yee's sample. Lawrence and 
Elvis (1981) have shown that the X-ray flux correlates with various optical 
and IR parameters for m/'st objects. However, NGC 6814 and NGC 4051 are under 
luminous In [Oil I], 3.5y and lOy flux relative their X-ray luminosities 
compared to the sample as a whole. 

We speculate that these correlations are consistent with the Idea that 
the X-ray emission began only a short time (' 100 yrs) ago for these 
objects. The probability of seeing one object In about 40 In the process of 
turning on depends greatly on the totally unknown X-ray lifetime of these 
objects. Consider the following possibility. After the X-ray source turns 
on. It will take some time x = R/c to Illuminate the entire broad line 
region. This will be the turn on time. The clouds quickly heat up and 
accelerate to v/c ' 1/50 and proceed to escape from the nuclear region. The 
escape time will be 50x. After the clouds escape both the X-ray source and 
the broad lines turn off. Thus one expects to see one object In ^ 50 In the 
process of turning on, which Is consistent with our Interpretation. Of 
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course, if this Is true then it Implies that an active nucleus has gone 
through mar\y on-off cycles in the history of the galaxy. 

If one turns on an X-ray source in an originally "normal* nucleus then 
one expects first to detect the effects on material closest to the nucleus. 
Since the broad line region is often less than one pc in size, an external 
observer would see the entire region illuminated in only a few years. 

However, the narrow line region which is 100-1000 pcs across will only be 
partially illuminated for "new" objects. Thermal reradiation from dust 
exterior to the nucleus should turn on with the narrow line region. In this 
picture Seyfert II galaxies represent the class which is turning off. Thus a 
large, faction of the observed differences in Seyfert galaxies would represent 
different snapshots of a single process. However tempting such a model is we 
note it does have problems. Oust near the nucleus will be at the highest 
temperature and will be observed to turn on first. This is contrary to 
observations in that NGC 6814 does have lOy emission (Rieke 1978) indicating a 
cool temperature and NG£ 1068 (a Seyfert II) has an IR spectra requiring a 
dust temperature of 1000 K (Jones et al. 1977). 

B. Inverse Compton Reflection 

In the inverse Compton reflection model, low energy (soft) photons enter 
a region of energetic electrons. The soft photons inverse Compton scatter 
from the electrons to emerge from the cloud as X-rays. In this section, we 
will examine the inverse Compton reflection process (Lightman and Rybicki 
1980) and find a constraint on the temperature of the soft photon source. We 
will make no^ assumptions about the electron population i.e., the distribution 
can be either thermal or non- thermal. We will assume that the X-ray emitting 
plasma is quite small and stable. By stable we mean that any variability seen 
is due to variations in the soft photon source and not due to changes in the 
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plasma Itself. Tennant et al. (1981) using the results of Lightman and 
Rybicki (1979) pointed out that the lack of spectral change during the 
intensity variations observed in NGC 6814 is consistent with this 
interpretation. Thus the lack of rapid variability in most sources could be 
explained by the growth of another “stronger" but more constant source of soft 
photons. If the X-ray source is slowly heating up its environment as 
mentioned in the last section then thermal re-radiation could be the new, 
constant, strong source of soft photons. This is consistent with the observed 
deficient IR flux for the rapidly varying galaxies. 

We will now find the minimum temperature that a thermal source can have 
and still provide enough photons (for Comptonization) to generate the observed 
spectrum. If the soft photons are at a temperature of kT then the observed 
power law will extend from ' 3 kT out to energies determined by the 
temperature of the scattering cloud. For the case of NGC 6814 the total 
number of photons radiated in the Comptonized spectrum is given by 

N^ = 4.7 x 10’^ 4it d^/^° dE (4) 

which corresponds to 1.2 x 10^2 ^ (^) at a distance d of 16 Mpc. The 
number of black body photons Impinging onto the X-ray region Is given by 



4.5 X 10 


49 2 


_e!l ( 

100 sec ' 


kT.3 

TeY^ 


(5) 


where r^QQ is the radius of the X-ray emitting region In units of 100 It 
ScC. NgB Is maxImiOT when the dilution factor Is equal to unity. Since for 
Compton scattering photon number Is conserved, we set Nj^^Ngg and find that 
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kT > 4.5 r 


- 1/2 

100 


eV. 


If r^QQ < 1 then we find the thermal source must have kT > 4.5 eV (52,000 
K). Dust grains at this temperature would quickly vaporize. This leads us to 
conclude that thermal radiation from dust cannot provide enough photons to 
generate the observed spectrum for a small X-ray source. This problem Is 
serious for NGC 6814 where the "hot spots" which provide the soft photons must 
be small and few In number In order to account for the rapid, large amplitude 
variability seen. If these spots reside outside the X-ray region then the 
dilution factor must be very small and hence kT » 4.5 eV. 

C. Source Size 

To summarize. If the source of soft photons is thermal In nature and If 
some of the optical and/or IR emission comes from the X-ray plasma directly, 
then the X-ray cloud nwst have r^QQ » 1. Since there appears to be some 
correlation between 3.5m IR emission and X-ray flux for most active galaxies 
(Lawrence and Elvis 1981) and since a large source size Is consistent with the 
absence of rapid X-ray variability reported here, we will consider the 
possibility of large X-ray regions. In order for lOOOK blackbody photons to 
be the soft photon source we find that r^QO ^ ^ days). This size 

1s consistent with previous observations of variability (Marshall, Warwick and 
Pounds 1981). An X-ray plasma this large could be generated In one of two 
ways In the black hole accretion picture. Either the central object Is large, 
hence very massive, or else the X-rays come from a large region not directly 
related to the central object. 

If the X-rays come from a region only » 10 gravitational radii from the 
central object and If the general lack of variability on time scales of less 
than one day tells us anything about the size of the region then the central 
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object must have a mass of 10® Mq's or greater. Pounds (1979) speculated that 
all active galaxies contained a 10® central object. The Eddington limit 
for such an object Is 10^^ erg/sec which Is much greater than the luminosity 
of any object In our sample. Of course. If most active galaxies contain a 10® 

Mq object then one Is forced to explain why NGC 6814 Is so small. It Is 

possible that most active galaxies contain dead quasars (Lynden-Bell 1969) 
whereas objects like NGC 6814 might never have gone through a “quasar 11ke“ 
phase. 

If the X-ray plasma Is not confined to the central object then It could 
fill a large volume of space. In this picture electrons are heated via some 
unknown mechanism, perhaps In a small volume and then proceed to fill a volume 
of space about one light day across. This may occur In a two temperature disk 

model (Shapiro, Lightman and Eardley 1976). Since the virlal temperature of 

the protons Is much higher than the corresponding electron temperature an 
accretion disk can have two temperatures. If the electrons are not 
effectively cooled then Interactions with protons can heat them to 
temperatures greatly exceeding their virlal temperatures. Lightman (1982), 
considering the possibility of relativistic electron plasmas, listed several 
non-gravltatlonal ways to confine the electrons. Since It Is notoriously 
difficult to confine plasmas via electromagnetic forces. It Is entirely 
possible that the central object Is boiling off some matter. Thus the X-ray 
emission can come from two components, the region near the black hole and the 
extended volume. 

To further consider a two component model for the X-ray source let us 
assume that the total power produced by a galactic nucleus comes out In two 
forms, first the Immediate production of X-rays and some Initially 
unobserved power P. The latter could be In the form of relativistic electrons 
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as mentioned above or in y rays as in the Penrose photoproduction model 
(Lei ter 1980). P will slowly fill a reservoir with energetic electrons. When 
steacly state is reached, the luminosity of the reservoir will be <P>. 
Therefore, the total X-ray luminosity will be 

+ a <P> » + r <L^> 

where o is the fraction of the reservoir's luminosity which comes out as 
X-rays and r is a <P>/<L^>. Below we will assume r is constant and 
that a ' 1 . 

Since only will show rapid variability let us consider what happens 
to B when one adds a variable source to a source of constant intensity Iq. If 
we assume that the intensities of the variable component are uniformly 
distributed between 0 and Ll then b for the sum is given by 

a = o,/I » Al/i/3 (2 I + aIH. (6) 

1 o 

We will now assume that the reservoir does not vary on the time scales we have 
sampled and also that the X-rays produced near the central object are highly 
variable. If we set Iq = ^Lj> (the luminosity of the reservoir) and Al = 2 
<L^> (the range of luminosities for immediate X-ray production) then 

® * TTlT+rT 

We define F to be the fraction of the total X-ray flux coming from the 
compact-variable region which is given by l/{l+r) = /3 F. 

We can check the consistency of this model by assuming that F > 75% and 
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that the average flux 1$ 4 for NQC 6814. Thus* If the assun^tlon of a uniform 
distribution of Intensities Is correct, we would expect to see Uie source vary 
from ' 1 to ' 7. The HOC 6814 flux shown In Figure 4 varies from 2 to - 7 
cts/cm^-k sec which Is In rough agreement with our model. 

Since our upper limits typically lie In the range of 10% we find that F 
Is typically less than 17%. For NGC 4151 which has a 0 ^ 6% we find that 
F « 10%. Thus the constant component would have to grow by a factor of 10 
(relative to the variable component) In order for 0 to decline from the NGC 
6814 value to the level observed for NGC 4151. One possibility Is that NGC 
6814 has not filled Its reservoir and thus has not come to steady state. This 
Is unlikely since it Implies that when steady state Is reached the liMsInosIty 
will be 10 times what It Is now. As shown In Tennant et a1. (1981) such a 
high total luminosity. If It were variable, would clearly violate the Fabian 
and Rees relation (1979). However It Is unclear as to whether the Fabian and 
Rees relation applies In this case since we are talking about a steady state 
condition. If the X-rays from the variable component pass through the 
reservoir, electron scattering could reduce the amplitude of variability. If 
the electrons have the correct power law distribution, scattering will not 
greatly alter the spectrum. One could also argue that NGC 6814 does not have 
a reservoir for ’ome unknown reason or that the X-ray production 
efficiency a for ':he reservoir Is low. 

D. Shot Noise 

In the GuHber:, Ross and Fabian (1982) model a cloud Is heated via some 
unknown mechanism /ind then allowed to cool via Inverse Comptonizatlon. Since 
the cooling times are very short the observed X-ray spectrum Is a time 
average. Their calculated "averaged" spectrum Is In good agreement with 
observations of the X-ray spectra of active galaxies. However, if the X-ray 




( 
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source Is powered by discrete events, these events should give rise to low 
amplitude variability. To Illustrate we apply the shot noise model. In this 
model \ events per time Interval each rise to amplitude h and then 

9 

exponentially decey away with time scale t. Thus using I <■ xhr and ■ 

2 

Xh t/2 we construct 

e • 

Me set N * Xt which rrpresents the number of "simultaneous" events taking 
place. If we assume that t Is longer than the smallest bln size we used (see 
Sutherland, Welsskopf, and Kahn 1978 Appendix A, to see how Isa function 
of bln size) and that x Is greater than one shot per day, then we can 
use 6 from Table 1. For NGC 6814 N « 10 which says that at any one time on 
the avenge 10 clouds dominate most of the X-ray flux. A more typical value 
of 0 near 10% Implies that N 200, and for Cen A, N Is greater than 1000. 
Since Cen A, NGC 4151 and NGC 6814 all have roughly the same lianlnoslty, 
accounting for the lack of rapid variability In Cen A and NGC 4151 by 
Increasing the number of shots only works If the shots become much more 
numerous and as a result each event much less luminous. 

One way to reconcile the Gullbert, Ross and Fabian model with the lack of 
variability Is to assume that the heating and cooling are taking place In a 
continuous manner. This argues for a stable accretion flow. 

E. Stabilized Accretion Flow 

Here we point out that a lack of variability can be due to dynamic 
stability. It Is clear from our observations of Cen A that there are periods 
when this highly variable source Is quite stable. For example, if the source 
was In the process of doubling its flux in N days, we can apply equation 3 and 
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f1nd» In a 12 hour observation, that e • 1/(2 m/l 2). If the (Mittrvtd value 
of 0 near 2S was mainly due to a linear trend then the doubling time mutt be 7 
deys (or longer). Thus periods of stability tre (^served In addition to the 
variability previously reported In the literature. 

It Is possible that flares could be created whenever a cloud Is fed Into 
the :entra1 engine. If NQC 6814 Is young then there may be many clouds close 
to the central engine. As the envlroraaent heats up the clouds wcHild disperse 
and the accretion flow becomes a uniform stream of gas. Cowle, Ostrlker and 
Stark (1978) showed that objects radiating with L > 1% of the Eddington 
liaalnoslty Mill be highly variable and that less luminous objects will be more 
stable. However, their model falls to generate the rapid time scale observed 
for NGC 6814. They asswned that the central object wac in a homogenous gas 
and the shortest time scale they obtained was related to the sound travel time 
across the sonic radius. Faster time scales may require the accretion of 
clouds. Although the Idea of lumpy vs smooth flow Is Interesting, an honest 
treatment of the problem requires detailed work beyond the scope of this 
paper. 

VI . CONCLUSIONS 

HEAO'l observatlo".^ of 38 active galactic nuclei have showr that rapid 
X-ray variability Is a rare occurrence. Of the objects In our sample only 
6814 showed significant variability on time scales less than 10^ seconds. 

Three other objects, all low luminosity Seyfert galaxies: NGC 4151, IteC 3227 
and MCG 5-23-16, showed variations consistent with doubling times of •• 1 
day. Our upper limits on variability correspond to roughly cfj/I < lOt on t\m 
scales from 300 to 10^ sec and o^/I < 20% on timescales from 5 to 300 sec. 

For the brightest objects In our sample Oj/I < 5% and for the weakest ®j/I < 

20 %. 
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W« discuss tht constraints thest obstrvatlons placa on sctnarlos of X-r^y 
production In active galactic nuclei. If these sources consist of two 
components, one highly vcriable and one quiescent thtm the ratio of the 
quiescent to variable flux Is typically greater than 10. If the X-rqy 
emission Is due to 1 nverse*Compton scattering of a thermal (« 1000K) source of 
soft photons we find that the typical X-ray emitting region must be larger 
than 3 light days. 

The general lack of variability on short time scales Indicates that 1) 
the X-ray region Is typically a light day or more across or^2) the process 
producing X-ray radiation Is extremely stable. If the X-ray region Is large 
then we can Infer that either the central object has a mass of > 10^ 

(Implying that the luminosity Is typically < 10"^ of the Eddington Vmitl or 
the X-ray emitting plasma Is not gravitationally confined to a few 
Schwarzschlld radii. Our present data do not allow us to discriminate among 
these possibilities. We also speculate that the difference between the 
avera^ object and the two objects that show rapid variability, NGC 6814 
(Tennant et al. 1981) and NGC 4051 (Marshall et al. 1981) (but see Natllsky, 
Shrader and Tananbaum 1982) Is because these two are In a sense, "newer*' 
objects than the typical Seyfert. That Is, many of their properties can be 
explained if their X-ray emission began recently (t < 100 years). 
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The value of that one calculates Is related not only to Intrinsic 
source variability (real or upper limits) but also to many experimental 
parameters such as the nun^er of photons collected and the amount of 
background, therefore, we have developed a simple expression which will 
relate the observed value of to a measure of the "true" source variability. 

For most of the sources we observed the background contributed more 
counts than the source. This means that for the smallest bln size of 5.12 sec 
we had more than 50 total counts (source and background) per bln. This large 
rate implies that the statistical errors are likely to be Gaussian and thus we 
can estimate the total variance with 


^ W 

2 = J J 


(x.-<x>)^ 2 

j ^ j ^ 


(Al) 


2 2 
where Wj = l/oj and we have used the definition of x . Since our background 

rate is extremely stable we will assume the total variance has two components, 
one due to photon ncise and another due to source variability 
Substituting into equation Al and solving for 



1 ) 


(A2) 


The statistical errors are based on the total number of counts seen in that 
data bin, thus we estimate 





(A3) 


28 


which when substituted Into Equation A2 gives 

^ f V ' ^ (A4) 

j ^ 

For the case of no background and no data gaps one can set o ^ equal to 

p 

the mean source count rate and then equation A2 reduces to the expression used 
by Forman, Jones and Tananbaum (1976). For observations with a large 
background one sets equal to the total (source plus background) count 
rate. The resulting equation was used by Parsignault and Grind! ey (1978). 
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TAK.E 1 


n ) 

(2) 

(3) 

(4) 

(5) 

(0) 

(7) 

(8) 

(9) 

OBJECT 

RA-DEC 

DATE 

LENCTH 

FLUX 

L 

cr/l 

071 

078 

CEN A 

1322-428 

384.5 

4 

40.5 

0.57 

1 .4 

<0.7 

<1 .6 

CEN A 

1322-428 

567.7 

11 

36.6 

0.43 

2.8 

1 .8 

3.0 

CYC A 

1057+408 

488.0 

6 

8.5 

05 

5.0 

<3.3 

<0.8 

ESO M1-G5B 

1017-588 

466.4 

0 

3.4 

16 

<8.1 

<5.7 

<0.8 

H 1640-505 

1648-501 

620.7 

0 

4.2 

0.28 

11.3 

10.0 

1 .0 

1C 4320A 

1346-301 

573.8 

0 

5.1 

4.1 

0.7 

4.6 

1.1 

MCC 2-58-22 

2302-000 

523.2 

6 

4.3 

34 

<6.0 

<4.0 

<0.5 

MCC 5-23-16 

0045-307 

407.7 

6 

8.1 

1 .0 

7.0 

5.7 

2.0 

MCC 8-11-11 

0551+464 

643.7 

0 

3.5 

5.1 

0.1 

0.8 

1 .8 

MK 142 

1022+510 

600.4 

0 

1 .6 

11 

<21.7 

<13.6 

<0.7 

m 279 

1352+600 

670.4 

13 

1 .2 

3.0 

36.5 

16.7 

0.9 

m. 270 

1352+606 

607.4 

7 

1.2 

3.0 

<18.4 

<11.2 

<0.7 

MK 335 

0m3+1OO 

553.8 

13 

1.1 

2.4 

31 .2 

18.3 

0.0 

MK 421 

1067+385 

513.1 

0 

1 .8 

6.0 

20.0 

11.0 

0.8 

MK 421 

1057+385 

704.4 

7 

1.2 

4.0 

<27.6 

<17.4 

<0.6 

MK 464 

1353+3M 

555.8 

10 

1 .4 

13 

<14.4 

<8.0 

<0.0 

MK 501 

1652+308 

610.8 

6 

2.7 

10 

<10.4 

<0.0 

<0.7 

MK 500 

2041-100 

403.7 

6 

4.6 

20 

5.4 

3.7 

0.8 

MK 500 

2041-100 

653.5 

4 

5.3 

23 

<6.8 

<5.3 

<1.3 

MK sra 

2041-100 

657.4 

18 

5.4 

24 

<4.3 

5.3 

’ .3 

MK 5M 

2041-100 

660 3 

• 6 

4.6 

20 

7.3 

6.4 

1 .4 

MK 500 

2041-1M 

666.4 

12 

3.6 

16 

7.0 

6.0 

1.0 

MK 500 

0212-010 

564.7 

6 

1 .1 

2.8 

<10.8 

<11.5 

<0.6 

M 82 

0051+690 

474.1 

5 

2.5 

0.M7 

<7.4 

<3.0 

<0.5 

NCC 526A 

0121-353 

526.8 

6 

2.1 

2.4 

17 3 

<6.2 

<0.6 

NCC 031 

0225+31 1 

585.0 

6 

1 .5 

1.4 

<18.7 

<12.0 

<0.7 

rCC 2110 

K4O-075 

647.6 

5 

5.0 

0.87 

<5.3 

<3.3 

<0.8 

NCC 2002 

0043-141 

500.4 

3 

7.8 

1.4 

7.4 

<2.5 

<0.0 

NCC 3227 

1020+201 

602.0 

6 

1 .7 

0.064 

25.5 

20.3 

1.6 

NCC 3783 

1136-375 

370.6 

3 

3.5 

1 .0 

<8.4 

<5.5 

<0.9 

NCC 3783 

1136-375 

554.8 

13 

3.1 

0.80 

14.7 

8.7 

1.2 

NCC 4151 

1208+307 

340.5 

3 

10.1 

0.38 

4.7 

3.4 

1 .5 

NCC 4151 

1208+307 

524.1 

3 

10.3 

0.73 

5.8 

6.5 

5.7 

NCC 4151 

1208+307 

532.1 

4 

22.8 

0.80 

<4.6 

<4.2 

<2.3 

NCC 4151 

1208+307 

533.0 

3 

18.8 

0.71 

<2.8 

<2.1 

<1 .8 

NCC 5506 

1410-030 

574.0 

6 

5.0 

0.64 

0.0 

6. 1 

1.5 

NCC 5548 

1415+254 

560.9 

6 

5.7 

5.4 

<4.0 

<3.1 

<0.8 

ICC 5549 

1415+254 

736.0 

0 

4.3 

4.1 

<5.4 

<3.8 

<0.8 

NCC 6814 

1940-104 

483.0 

7 

4.1 

0.40 

43.6 

43.3 

8.4 

NCC 7213 

2206-474 

511 .4 

6 

4.2 

0.40 

<8.4 

<5.8 

<0.8 

NCC 7460 

2300+086 

547.0 

12 

1.8 

1 .7 

10.4 

16.0 

1.3 

NCC 7582 

2315-426 

518.3 

6 

3.1 

0.26 

0.2 

9.0 

1 .5 

PKS 0548-322 

0540-323 

638.0 

8 

2.6 

43 

<10. 1 

<8.5 

0.7 

PKS 2151-304 

2151-306 

677.7 

5 

8.1 

810 

<5.2 

5.1 

1 .6 

2A 1210-^305 

1210+305 

516.2 

6 

3.0 

180 

<8.5 

<5.4 

<0.8 

3C 111 

0415+370 

611.2 

0 

3.0 

20 

<8.3 

<6.4 

<1 .2 

3C 120 

0430+052 

615.6 

0 

3.1 

12 

<7.2 

<4.7 

<0.7 

3C 273 

1226+023 

533.0 

6 

0.4 

800 

4.8 

3.3 

1.5 

3C 273 

1226+023 

546.4 

6 

8.5 

720 

<3.3 

<2.3 

<0.0 

3C 273 

1226+023 

540.3 

8 

0.0 

7M 

<3.4 

<2.4 

<1 .0 

3C 273 

1226+023 

551 .2 

10 

8.7 

740 

<3.3 

<2.4 

<1 .0 

3C 382 

1833+327 

601 .5 

0 

3.2 

38 

12.2 

<6.1 

<0.0 

3C 300.3 

1845+707 

723.2 

10 

2.2 

25 

<14.1 

<0.4 

<0.8 


NOTES i 

(1) AlphabsTlc by source nams 

(2) RA (hours and mins) Osc (degress and tenths) 

(3) Beginning of obeervatlont day of year 1077 
(1) Length of observation in hours 

(5) Cts/c»2-||««c 

(6) 2-20 Kev luminosity (H-7K km/sec/Hpc) asscmlng no absorption » unlts-10^* ergs/ssc 

(7) 328s bln size 

(8) 86n bln size 


(0) 86m bln size 
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FIGURE CAPTIONS 


Figure 1 - Shows e, a measure of Intrinsic source varlablll^ (eq. 2), vs 
source flux. The data was collected Into 'o 5 min bins In order to 
construct this figure. The line Is an estimate of the error Introduced 
by source confusion (see text). The solid dots are positive detections 
of source variance while the bars are upper limits (see Table 1). 

Figure 2 - Same as Figure 1 except now the bins are about one spacecraft 
orbit long. 

Figure 3 top - A histogram of Oj divided by the background counting rate for 
our observations of active galaxies. The hatched regions represent 
positive detections of variability whereas the white area above the 
hatched region represent the 90% confidence upper limit for the 
non-variable objects. 

bottom - Same as top panel except In this case no HEAO-1 sources were 
In the field of view. 

Figure 4 - A sample of 9 of our X-ray light curves. The top 6 panels 
represent the more variable objects In our sample whereas the bottom 3 
represent typical light curves. In all cases a 20 min bln size was used 
to construct the light curves. The first number Is Oj/I and the nwnber 
In parentheses Is Oj/B. Both Oj/I and Oj/B are constructed for the bln 
size given In the plot. The best fitting linear trend Is Indicated by a 
solid line for MCG 5-23-16 and NGC 3227. 

Figure 5 - The most variable observation of M6C 4151 has been fitted to 
three models: (top) a straight line through all the points* (middle) a 
step function, and (bottom) a straight line through the last 5 points. 
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